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ABSTRACT 

We report a systematic and statistically significant offset between the optical (g—z or B-I) colors of seven 
massive elliptical galaxies and the mean colors of their associated massive metal-rich globular clusters (GCs) in 
the sense that the parent galaxies are redder by ~ 0.12-0.20 mag at a given galactocentric distance. However, 
spectroscopic indices in the blue indicate that the luminosity-weighted ages and metallicities of such galaxies 
are equal to that of their averaged massive metal-rich GCs at a given galactocentric distance, to within small 
uncertainties. The observed color differences between the red GC systems and their parent galaxies cannot be 
explained by the presence of multiple stellar generations in massive metal-rich GCs, as the impact of the latter 
to the populations' integrated g—z or B-I colors is found to be negligible. However, we show that this paradox 
can be explained if the stellar initial mass function (IMF) in these massive elliptical galaxies was significantly 
steeper at subsolar masses than canonical IMFs derived from star counts in the solar neighborhood, with the 
GC colors having become bluer due to dynamical evolution, causing a significant flattening of the stellar MF 
of the average surviving GC. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: stellar content — galaxies: formation — 
galaxies: star clusters: general — globular clusters: general 



1. INTRODUCTION 

The shape of the stellar initial mass function (IMF) is a fun- 
damental property for studies of star and galaxy formation. It 
constitutes a crucial assumption when deriving physical pa- 
rameters from observations. However, its origin remains rela- 
tively poorly understood. A particularly important question is 
whether or not the IMF is "universal" among different envi- 
ronments, since galactic mass-to-light (A4/L) ratios are very 
sensitive to variations in the shape of the IMF at subsolar 
masses. Several observational studies find that the IMF in 
various environments within our Galaxy and the Magellanic 
Clouds seems remarkably uniform down to detection limits 
of - 0.1 - 1M Q (Kroupa 2001; Chabrier 2003; Bastian et al. 
2010, and references therein). However, recent studies indi- 
cate that the IMF may be different in more extreme environ- 
ments: strong gravity-sensitive features in near-IR spectra of 
luminous elliptical galaxies by van Dokkum & Conroy (2010, 
201 1) favor a bottom-heavy IMF (a ~ -3 in dN/dM oc M a ). 
If confirmed, this result would have widespread implications 
on studies of stellar populations and galaxy evolution. For 
example, such a bottom-heavy IMF would increase the M. /L 
ratio in the K band by a factor ~ 3 relative to a "standard" 
Kroupa (2001) IMF, and hence invalidate the widespread as- 
sumption that M /L ratios in the near-IR are highly indepen- 
dent of galaxy type or mass. 

In this paper, we describe an effort to find independent evi- 
dence to confirm or deny the presence of a bottom-heavy IMF 
in nearby giant elliptical galaxies. We use observed optical 
colors and spectral line indices of metal-rich globular clusters 
(GCs) and their parent galaxies in conjunction with dynamical 
evolution modeling of GCs. 

Infrared studies of star formation within molecular clouds 
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have shown that stars form in clusters or unbound associa- 
tions with initial masses A^ c i,o m the range 10 2 - 10 8 M Q 
(e.g., Lada & Lada 2003; Portegies Zwart et al. 2010; Kmijs- 
sen 2012, and references therein). While most star clusters 
with A^ci.o S 10 4 Mq are thought to disperse into the field 
population of galaxies within a few Gyr, the surviving mas- 
sive GCs constitute luminous compact sources that can be ob- 
served out to distances of several tens of megaparsecs. Fur- 
thermore, star clusters represent very good approximation s of 
a "simple stellar population" (hereafter SSP), i.e., a coeval 
population of stars with a single metallicity, whereas the dif- 
fuse light of galaxies is typically composed of a mixture of 
different populations. Thus, star clusters represent invaluable 
probes of the star formation rate (SFR) and chemical enrich- 
ment occurring during the main star formation epochs within 
a galaxy's assembly history. 

Deep imaging studies with the Hubble Space Telescope 
(HST) revealed that giant elliptical galaxies typically contain 
rich GC systems with bimodal color distributions (e.g., Kundu 
& Whitmore 2001; Peng et al. 2006). Follow-up spectroscopy 
showed that both "blue" and "red" GC subpopulations are 
nearly universally old, with ages >8 Gyr (e.g., Cohen et al. 
2003; Puzia et al. 2005; Brodie et al. 2012). This implies that 
the color bimodality is mainly due to differences in metal- 
licity. The colors and spatial distributions of the blue GCs 
are usually consistent with those of metal-poor halo GCs in 
our Galaxy and M31, while red GCs have colors and spatial 
distributions that are similar to those of the "bulge" light of 
their host galaxies (e.g., Geisler et al. 1996; Rhode & Zepf 
2001; Bassino et al. 2006; Brodie & Strader 2006; Peng et al. 
2006; Goudfrooij et al. 2007). Thus, the red metal-rich GCs 
are commonly considered to be physically associated with the 
stellar body of giant elliptical galaxies. This represents an as- 
sumption for the remainder of this paper. 

2. GALAXY SAMPLE AND GC SELECTION 

Our main galaxy selection criterion is the presence of a very 
clear bimodal optical color distribution of its GCs as derived 
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from high-quality, high spatial-resolution (i.e., HST) imaging 
that allows GC detection as close to the galaxy centers as pos- 
sible. With this in mind, most galaxies in our sample were 
drawn from the ACS Virgo Cluster Survey (ACSVCS; Cote 
et al. 2004), a homogeneous survey of 100 early-type galax- 
ies and their GCs in the Virgo galaxy cluster, using the F475 W 
and F850LP filters (hereafter g and z, respectively). Specifi- 
cally, we use galaxy photometry from Ferrarese et al. (2006) 
supplemented by GC data from Jordan et al. (2009). Galaxies 
were selected based on the properties of GC color distribu- 
tions as modeled by Peng et al. (2006) using the Kaye's Mix- 
ture Model (KMM; McLachlan & Basford 1988; Ashman et 
al. 1994), fitting two Gaussians to the data. Two criteria were 
used for selection: (i) the p-value of the KMM fit needs to 
be < 0.01, and (ii) the color distribution of the GCs needs 
to clearly show two distinct peaks as judged by eye in Fig. 2 
of Peng et al. (2006). Finally, we exclude NGC 4365 due to 
conflicting results on the presence of intermediate-age GCs 
among its red GCs in the literature 3 . This selection yields the 
high-luminosity galaxies NGC 4472, NGC 4486, NGC 4649, 
NGC 4552, NGC 4621, and NGC 4473 (see Table 1). 

In the context of the current study, it is important to include 
a comparison between metallicity measurements of red GCs 
from colors versus from spectral line indices (see Section 3 
below). As adequate line index data for several red GCs is 
currently only available for one galaxy in the ACSVCS sam- 
ple, we add NGC 1407 to our galaxy sample given that high- 
quality HST/ACS imaging (in filters F435 W (hereafter B) and 
F814W (I)) as well as optical spectroscopy of both red GCs 
and the diffuse light are available (Forbes et al. 2006; Cenarro 
et al. 2007; Spolaor et al. 2008a,b). 

The selection of red GCs in these galaxies was done by se- 
lecting GCs redder than the color for which the KMM proba- 
bility of a GC being member of the blue versus the red peak 
is equal. In addition, we require GCs to be massive enough 
to render the effect of stochastic fluctuations of the number of 
red giant branch (RGB) and asymptotic giant branch (AGB) 
stars on the integrated colors smaller than the typical photo- 
metric error of a(g-z) = 0.04 mag. Using the methodology 
of Cervino & Luridiana (2004) for an age of 12 Gyr and solar 
metallicity 4 , this requirement translates in a minimum cluster 
mass Aid — 2.5 x 10 5 M Q . For galaxies in the Virgo cluster, 
this mass limit is equivalent to selecting GCs with zab < 21 .91 
mag under the assumptions of (i) M c \/L z =1.5 M©/L Z!0 for 
Z = Z© (this value varies only slightly with metallicity; Jordan 
et al. 2007) and (ii) m-M = 31.09 mag for the Virgo clus- 
ter (Jordan et al. 2009). The equivalent magnitude limit for 
the GCs in NGC 1407 is / < 23.28 mag, using M C \IU = 2.1 
Mq I Li Q for Z = Z© (Maraston 2005) and m-M = 3 1 .60 mag 
(Forbes et al. 2006). 

3. RADIAL COLOR GRADIENTS OF GALAXIES AND THEIR RED 
GC SYSTEMS 

In this section we compare the colors of the red GCs (se- 
lected as mentioned above) with those of the diffuse light of 
their parent galaxies as a function of projected galactocentric 
radius (hereafter R gd i). All observed colors for the ACSVCS 
galaxies were dereddened using the Galactic foreground red- 
dening values listed in Jordan et al. (2009). The colors of GCs 

3 See Puzia et al. (2002), Larsen et al. (2003), and Kundu et al. (2005) 
versus Brodie et al. (2005) and Chies-Santos et al. (201 1). 

4 Fluctuation magnitudes for the g and z passbands are estimated by inter- 
polating between those of B and V and 1q and J, respectively, using the values 
in Table 5B of Worthey (1994). 



in NGC 1407 from Forbes et al. (2006) were already corrected 
for Galactic foreground reddening. For each galaxy we cal- 
culate a mean color of the red GC population along with its 
radial gradient by means of a weighted linear least-squares fit 
between the GCs colors col (i.e., g—z or B-I) and the loga- 
rithm of R gd \ in arcsec: 

col (R g . dl ) = col + G cot log(/? gal ) ( 1 ) 

where colo is the g—z or B-I color at log (R ga \) = and G co i is 
the color gradient per dex in radius, A( , co/)/Alog(/? ga i). For 
each GC, uncertainties associated with its measured color and 
its membership of the red population are taken into account 
in the fit by adding the following two parameters in quadra- 
ture: (7) the photometric error of the GCs color, and (ii) the 
inverse probability of the GC being a member of the red popu- 
lation (p re d)» as calculated by the KMM algorithm mentioned 
above, expressed in magnitudes (i.e., log(l//? rec i)/log(2.5)). 
Radial color profiles for the diffuse light of the parent galax- 
ies were determined using the same functional form as Eq. (1). 
For NGC 1407, B-I data of the diffuse light were taken from 
the surface photometry of Spolaor et al. (2008a), which was 
derived from the same HST/ACS data as the GC photometry 
of Forbes et al. (2006). ACSVCS galaxy color gradients are 
taken from Liu et al. (201 1), who used the surface photometry 
of Ferrarese et al. (2006). Observed ACSVCS galaxy colors 
at Rga\ = 1 arcsec were taken from Ferrarese et al. (2006) and 
listed in Table 1 . Panels (a) - (c) and (e) - (g) of Figure 1 show 
g—z versus R gd \ for the 6 ACSVCS galaxies and their massive 
red GCs. Panel (a) of Figure 2 shows the equivalent for NGC 
1407 and its red GCs. Note that while the radial gradients 
of the galaxies' colors are similar to those of the "average" 
colors of their red GC systems, there is a systematic offset be- 
tween the two sets of radial color profiles in the sense that the 
average color of red GCs is ~ 0. 12 - 0.20 mag bluer in g—z or 
B-I than that of their parent galaxies. This color difference 
was already noted by Peng et al. (2006), but they did not pur- 
sue an investigation of its possible cause(s). The significance 
of this color difference is 7 a to 22 a depending on the galaxy, 
where a is the mean error of the fit of eq. (1) to the colors of 
red GCs in the sample galaxies (see Table 1). 



TABLE 1 

Colors of sample galaxies and their massive red GCs. 



NGC 




(g-z)o.gal 


fe-z)o 


G s-z 


rms 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


4472 


-22.90 


1.586 ±0.005 


1.425 


-0.021 


0.010 


4486 


-22.66 


1.598 ±0.005 


1.413 


+0.000 


0.007 


4649 


-22.41 


1.592 ±0.005 


1.449 


-0.022 


0.010 


4552 


-21.36 


1.558 ±0.005 


1.372 


+0.000 


0.014 


4621 


-21.18 


1.596 ± 0.005 


1.365 


-0.023 


0.013 


4473 


-20.70 


1.563 ± 0.005 


1.358 


-0.005 


0.017 


NGC 




(B-/)o,gai 


(B-Do 


Gb-i 


rms 


1407 


-21.86 


2.313 ±0.005 


2.103 


-0.004 


0.007 



NOTE. — Column (1): NGC number of galaxy. Column (2): absolute V 
magnitude of galaxy. Column (3:) colo of galaxy. Column (4): colo of red 
GCs. Column (5): color gradient G co i of red GCs in mag/dex. Column (6): 
mean rms error of fit of Eq. (1) to colors of red GCs (in mag). 



Under the assumption that surviving GCs represent the 
high-mass end of the star formation process that also popu- 
lated the field (Elmegreen & Efremov 1997), this color dif- 
ference is difficult to understand in terms of a metallicity dif- 
ference. While we deem it likely that the spread of colors 
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FIG. 1 . — Radial color gradients of giant elliptical galaxies from the ACSVCS survey and their systems of massive red GCs. Panels (a) - (c) and (e) - (g) show 
dereddened (#-z)ab versus projected galactocentric radius for the ACSVCS galaxies mentioned at the top left of the panel in question. Black solid lines represent 
linear least-square fits to the galaxy colors, red circles with black error bars represent individual red GCs, and the red dashed lines depict linear least-square fits 
to the red GC systems as a whole. The dashed lines in panel (d) plot integrated (g-z)AB versus mass function (MF) slope olmf for Marigo et al. (2008) isochrones 
for an age of 12 Gyr and the two [Z/H] values shown in the legend. The dotted lines in panel (d) indicate the colors for a Kroupa (2001) MF. Panel (h) is similar 
to panel (d), but now the upper (black) lines are drawn for [Z/H] = -0.2. The grey shaded region crossing panels (a) through (d) and (e) through (h) illustrate the 
range of MF slopes reproducing the galaxy colors for the upper isochrone lines drawn in panels (d) and (h), respectively. 
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Fig. 2. — Similar to Figure 1, but now for NGC 1407 and its red GCs. 
Panel (a) shows (B—I)q versus projected galactocentric radius. The dashed 
lines in panel (b) plot integrated (B—I)o versus MF slope olmf for Marigo et 
al. (2008) isochrones for an age of 12 Gyr and the two [Z/H] values shown in 
the legend. Symbols, lines, and shading as in Figure 1. 



exhibited by the red GC systems in Figure 1 is due (at least in 
part) to a range of GC metallicities (and/or ages), one would 
expect that the same range of metallicities (and/or ages) be 
present among the field stars and hence that the average col- 
ors would be similar. 

This argument is corroborated by spectroscopic Lick in- 
dex data for red GCs and the diffuse light in the two gi- 
ant elliptical galaxies in our sample (or in the literature at 
large) for which we found such data available in the litera- 
ture for both a significant number of red GCs and the dif- 
fuse light within a range in R i3 \ sampled by both compo- 
nents. As to the GC spectra, Cohen et al. (2003) obtained 
Keck/LRIS spectra of 47 GCs in NGC 4472, while Cenarro 
et al. (2007) used Keck/LRIS to obtain spectra for 20 GCs in 
NGC 1407. In both studies, the resulting GCs were split ap- 
proximately evenly between blue and red GCs. To evaluate 
spectroscopic ages and metallicities, we use the indices H/3 

and [MgFe]' = (Mgb x (0.72 Fe 52 70 + 0.28 Fe 53 35)) 1/2 , respec- 
tively, and compare them to predictions of the SSP models 
of Thomas et al. (2003). The [MgFe]' metallicity index was 
chosen because Thomas et al. (2003) found it to be indepen- 



dent of variations in [a/Fe] ratio. This is relevant since giant 
elliptical galaxies and their metal-rich GCs are known to ex- 
hibit supersolar a/Fe ratios (e.g., Trager et al. 2000; Puzia et 
al. 2006). Figure 3 shows H/3 versus [MgFe]' for the GCs 
in NGC 4472 from Cohen et al. (2003), subdivided into blue 
and red GCs. As the uncertainties are significant for the in- 
dividual GCs, we also indicate weighted average indices for 
the blue and red GCs (see large blue and red symbols in Fig- 
ure 3). Note that [Z/H] ~ +0.30 for the average red GC in 
Figure 3, which corresponds to [Z/H] = -0.07 on the Zinn & 
West (1984) scale (Cohen et al. 2003). For comparison, we 
overplot H/3 and [MgFe]' for the diffuse light of NGC 4472 
(Davies et al. 1993; Fisher et al. 1995) at # ga i = 50". Note that 
the median R gSL \ of the red GCs targeted by Cohen et al. (2003) 
is 59"1. Similarly, Figure 4 shows H/3 versus [MgFe]' for the 
blue and red GCs in NGC 1407 from Cenarro et al. (2007) and 
the diffuse light of NGC 1407 at fl ga i = 50" from Spolaor et 
al. (2008b). For comparison, the median R gd i of the red GCs 
targeted by Cenarro et al. (2007) is 50. "3. We conclude that 
the average red GC and the diffuse light have equal values of 
H/3 and [MgFe]' to within la for both NGC 1407 and NGC 
4472, the two galaxies for which we could find spectroscopic 
ages and metallicities for both red GCs and the diffuse light, 
at comparable galactocentric distances, in the literature. This 
indicates that the mean ages and metallicities of the red GCs 
and the underlying diffuse light of giant elliptical galaxies are 
indeed the same within the uncertainties. 

We thus arrive at a situation where the mean spectroscopic 
ages and metallicities of the diffuse light of giant elliptical 
galaxies and their red GCs seem consistent with one another 
while the g-z or B-I colors of the diffuse light of the galaxies 
are systematically redder than those of their average red GC . 
We evaluate three possible solutions to this paradox below. 

4. POSSIBLE CAUSES OF THE COLOR OFFSET BETWEEN RED GC 
SYSTEMS AND THEIR PARENT GALAXIES 

4.1. A Mismatch between Photometric and Spectroscopic 
Samples? 
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FIG . 3 .— H/3 versus [MgFe]' for GCs in NGC 4472 from the data of Cohen 
et al. (2003). Small blue and red circles indicate GCs designated by Geisler et 
al. (1996) as "blue" and "red", respectively. Large filled blue and red squares 
represent weighted average values of H/3 and [MgFe]' of the blue and red 
GCs, respectively. Large open red circles indicate red GCs with HST/ACS 
photometry (see Figure 5a and Sect. 4.1). Error bars at the top right indi- 
cate typical uncertainties for individual GCs. SSP models of Thomas et al. 
(2003) are overplotted for [o/Fe] = +0.3. Dashed lines indicate [Z/H] values 
of -2.25, -1.35, -0.33, 0.00, 0.35, and 0.67 dex. Solid lines indicate ages of 
1, 2, 3, 5, 10, and 15 Gyr. For comparison, the black triangle represents H/3 
and [MgFe]' for the diffuse light of NGC 4472 at R gal = 50". See Sect. 3. 
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FIG. 4. — Similar to Figure 3, but now for GCs in NGC 1407 from the data 
of Cenarro et al. (2007) and for the diffuse light of NGC 1407 at i? gal = 50" 
from the data of Spolaor et al. (2008b). Large open red circles are not shown 
in this case since all GCs shown have HST/ACS photometry (cf. Figure 5b). 

As a first possible solution to the paradox described in the 
previous Section, we consider the hypothesis that the red GCs 
in NGC 4472 and NGC 1407 whose spectroscopic Lick index 
data were shown in Figures 3 and 4, respectively, happen to be 
redder and more metal-rich on average than the respective full 
photometric samples of red GCs. This hypothesis is rejected 
by the data, as illustrated in Figure 5 in which open circles 
highlight the red GCs that have counterparts in Figures 3 and 



4. Dotted black lines indicate the mean g—z or B— I colors 
of the GCs with Lick index data, which are consistent with 
the mean colors of the full red GC systems to within 1 a. Re- 
calling that the significance of the observed color differences 
between the metal-rich GC systems and their host galaxies is 
in the range 7 - 22 a, we conclude that the red GCs with Lick 
index data are not misrepresenting the full systems of red GCs 
in terms of their g—z or B—I colors. 

— i — i i ■ ■ ■ ■ 1 1 — i — i i ■ ■ ■ ■ 1 1 




10 1 10 a 
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FIG. 5. — Panels (a) and (b) are near-copies of Figures la and 2a, respec- 
tively, except that error bars are now shown in grey. Large open circles in 
panels (a) and (b) identify the red GCs in the HST data for which Lick in- 
dex data were shown in Figures 3 and 4, respectively. The black dotted lines 
indicate the mean colors of the two sets of red GCs shown with large open 
circles. Note that the black dotted lines almost coincide with the red dashed 
lines, which depict linear fits to Eq. 1 for the full red GC systems. 



4.2. An Effect of the Presence of Multiple Stellar Populations 
in Massive GCs? 

Several detailed photometric and spectroscopic investiga- 
tions over the last decade have established that massive Galac- 
tic GCs typically host multiple, approximately coeval, stel- 
lar populations (see the recent review by Gratton et al. 2012, 
and references therein). The recognition of this fact came 
mainly with the common presence of anticorrelations between 
[Na/Fe] and [O/Fe] among individual stars (often dubbed 
"Na-O anticorrelations"; see, e.g., Carretta et al. 2010), which 
are seen within virtually all Galactic GCs studied using multi- 
object spectroscopy with 8-m-class telescopes to date. The 
leading theory on the cause of these abundance variations is 
the presence of a second generation of stars formed out of 
enriched material lost (with low outflow velocity) by "pol- 
luters" (massive stars and/or intermediate-mass AGB stars) 
of the first generation. High-temperature proton-capture nu- 
cleosynthesis in the atmospheres of the "polluters" increased 
the abundances of He, N, and Na relative to Fe (e.g., Renzini 
2008). 

Note that this theory only predicts observable star-to-star 
abundance variations among light elements (up to 13 A1), 
which is consistent with what is seen among almost all Galac- 
tic GCs. Variations in heavy (e.g., Fe-peak) element abun- 
dances have only been found within two very massive metal- 
poor GCs: wCen (e.g. Lee et al. 2005; Johnson & Pila- 
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chowski 2010) and M22 (Marino et al. 2009). It is thought 
that such GCs represent surviving nuclei of accreted dwarf 
galaxies whose escape velocities were high enough to retain 
Fe-peak elements from SN explosions (e.g., Lee et al. 2007; 
Georgiev et al. 2009). However, we deem this accretion sce- 
nario rather unlikely for the formation history of red GCs 
given their high metallicities. We therefore restrict the fol- 
lowing discussion to GCs without a spread in iron abundance. 

The effect of light-element abundance variations within 
GCs to photometry seems to be more subtle than that to 
spectroscopic line indices. While photometric splitting of 
RGB, main sequence (MS), and/or subgiant branch (SGB) 
sequences have been observed in color-magnitude diagrams 
(CMDs) of several Galactic GCs, such splits only show up in 
CMDs that involve filters shortward of ~ 4000 A, while they 
disappear when using visual passbands such as B, V, and/or / 
(see, e.g., Marino et al. 2008; Kravtsov et al. 201 1; Milone et 
al. 2012). 

To evaluate the effect of the presence of multiple stellar 
populations in massive GCs to their colors, we determine col- 
ors of various stellar types from synthetic spectra with chemi- 
cal compositions typical of first- and second-generation stars 
found in massive GCs, and we use the resulting stellar spec- 
tra to calculate integrated colors of the full stellar popula- 
tion. This procedure is described in detail in Appendix A. 
As shown there, the expected impact of the presence of two 
stellar generations to the integrated B-I and g— z colors of 
metal-rich GCs is offsets of -0.021 and -0.027 mag, respec- 
tively, relative to a true SSP. As this is only a small fraction 
of the observed color offsets between red GC systems and the 
diffuse light of giant elliptical galaxies, it seems that an addi- 
tional effect is needed to produce that color offset. 

4.3. A Steep Stellar Mass Function Slope in the Field Star 
Component? 

At a given age and chemical composition, the only funda- 
mental property of a simple stellar population that can affect 
its integrated color (or spectral line index) significantly is the 
shape of its stellar MF. Figures Id, lh, and 2b illustrate the ef- 
fect of changing the slope a of the stellar MF to the integrated 
g— z or B-I color of a stellar population. The dashed curves 
in these panels were determined from Padova isochrones 
(Marigo et al. 2008) of age 12 Gyr and for the [Z/H] values 
shown in the figure. After rebinning the isochrone tables to a 
uniform bin size in the stellar mass M.* using linear interpo- 
lation, weighted luminosities of the population were derived 
for the relevant passbands by weighting individual stellar lu- 
minosities by a factor (Ai+/ ' M. m ?a) a where A4 mdx is the max- 
imum stellar mass reached in the isochrone table. For refer- 
ence, the dotted lines show the expected colors for a "stan- 
dard" Kroupa (200 1 ) IMF. 

A comparison of panels (a) - (c) with panel (d) of Figure 1 
reveals several items of interest. First, the "average" colors of 
the red GC systems are consistent with the prediction of a SSP 
model with age =12 Gyr, [Z/H] = 0.0, and a Kroupa MF. This 
is also consistent with the spectroscopic age and [Z/H] found 
for the average red GC found above in Section 3. Second, the 
redder colors of the parent galaxies with respect to the average 
colors of the red GC systems can be explained if the MF of the 
field star component is "bottom-heavy" (-3.0 < a < -2.4 for 
[Z/H] = 0.0) relative to that of the red GC system, especially 
in the inner regions. This assumes that the field stars and 
the massive red GCs share the same distributions of age and 



metallicity, for which we showed spectroscopic evidence in 
Figs. 3 and 4. In this regard it is important to note that the de- 
pendence of the blue-visual spectral indices H/3 and [MgFe]' 
to changes in the MF slope in the range indicated by the color 
differences between red GCs and their parent galaxies is neg- 
ligible. This is illustrated in Figure 6 where we compare the 
values of H^ and [MgFe]' for NGC 4472 and NGC 1407 at 
R ga i = 50" (cf. Figs. 3 and 4) with predictions of the Vazdekis 
et al. (2010) SSP models which were calculated for a variety 
of IMF types and slopes 5 . We plot model predictions for (i) a 
Kroupa IMF, (ii) an IMF with a Salpeter slope (a = -2.3) for 
stellar masses M > 0.6 M and a = -2.8 for M < 0.6 M Q , 
and (Hi) an IMF with a = -2.8 for all stellar masses. As Fig- 
ure 6 shows, the three sets of model predictions are consistent 
with one another to within the measurement uncertainties in 
the relevant region of parameter space. 



~> I ' ' ' I ' 1 1 I 1 

6 - Vazdekis et al. (2010) models - 



-2.32 dex 




+ 0.22 



I I I I I I I I I I I I I 

2 4 

[MgFe]' (A) 

FIG. 6. — H/3 versus [MgFe]' for SSP model predictions of Vazdekis et al. 
(2010) using different IMFs (see legend on the right). Black lines: Kroupa 
(2001) IMF. Red lines: a = -2.8 for M < 0.6 Mq and Salpeter IMF for 
M > 0.6 Mq. Blue lines: a = -2.8 for all stellar masses. Dashed lines 
indicate [Z/H] values of -2.32, -1.71, -1.31, -0.71, -0.40, 0.00, and 0.22 
dex. Solid lines indicate ages of 1, 2, 3, 5, 10, and 15 Gyr. For comparison, 
the grey and open triangles represent H/3 and [MgFe]' at Rg^ = 50" for the 
diffuse light of NGC 4472 and NGC 1407, respectively. See the discussion 
in Sect. 4.3. 

Similar conclusions can be drawn when comparing pan- 
els (e)-(g) with panel (h) of Figure 1 and panel (a) with 
panel (b) of Figure 2, except that the average colors of the 
red GC systems of the less luminous galaxies depicted there 
(see Table 1) are better matched by a slightly lower metal- 
licity ([Z/H] ~ -0.2). This is consistent with the color- 
magnitude relation among elliptical galaxies when interpret- 
ing color changes in terms of metallicity changes: the slope in 
the V—K vs. V relation among E galaxies in the Virgo galaxy 
cluster is -0.079 ±0.007 (Bower et al. 1992), which translates 
to a [Z/H] vs. galaxy magnitude relation of -0.093 ± 0.008 
dex/mag according to the Marigo et al. (2008) SSP models at 
an age of 12 Gyr. 

5 Note that Vazdekis et al. (2010) defined IMF slopes fj, as dN/dM oc 
Hence, our a = -(p,+ \ ). 
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A possible mechanism for producing the steep MFs in giant 
elliptical galaxies (relative to those of their massive red GCs) 
indicated by our results is discussed in the next Section. 

5. DYNAMICAL EVOLUTION EFFECTS ON STELLAR MASS 
FUNCTIONS OF GLOBULAR CLUSTERS 

The gradual disruption of star clusters over time affects the 
shape of their stellar MF, because the escape probability of 
stars from their parent cluster increases with decreasing s tellar 
mass (e.g., Baumgardt & Makino 2003). This effect has been 
used to explain the observed stellar MFs in ancient Galactic 
GCs, some of which are flatter than canonical IMFs (e.g., de 
Marchi et al. 2007; Paust et al. 2010). Such flat MFs in GCs 
can explain their low observed mass-to-light ratios when com- 
pared with predictions of SSP models that use canonical IMFs 
(e.g., Kruijssen & Mieske 2009). The flattening of the MFs in 
Galactic GCs is greatest for GCs that are closest to dissolution 
(Baumgardt et al. 2008; Kruijssen 2009). 

Here, we explore the hypothesis that the IMF in the (high- 
metallicity, and likely high-density) environments that pre- 
vailed during the main star formation events that formed the 
stars in the bulge and the red GCs of giant elliptical galax- 
ies was relatively steep (-3.0 < a < -2.4, as indicated by 
the analysis in the previous Section), i.e., steeper than the 
canonical IMF seen in the solar neighborhood, and that dy- 
namical evolution of the star clusters formed during those 
events caused the stellar MF of the "average" surviving mas- 
sive star cluster to be similar to a canonical IMF. To this end, 
we use the GC evolution model of Kruijssen (2009, hereafter 
K09) which incorporates the effects of stellar evolution, stel- 
lar remnant retention, dynamical evolution in a tidal field, and 
mass segregation. For the purposes of the current study, we 
select K09 models that ( i) feature solar metallicity and "de- 
fault" kick velocities of stellar remnants (white dwarfs, neu- 
tron stars, and black holes), and ( ii) produce GCs with masses 
M d >2.5 x 10 5 M Q at an age of 12 Gyr, i.e., the GCs shown 
in Figures 1 and 2 (cf. Section 2). We further consider King 
(1966) profiles with Wo values of 5 and 7 and environmen- 
tally dependent cluster dissolution time scales fo 6 of 0.3, 0.6, 
1 .0, and 3.0 Myr. For reference, f o = 1.3 Myr yields a reason- 
able fit to the globular cluster mass function of all (surviving) 
Galactic GCs (Kruijssen & Portegies Zwart 2009). 

Note however that GC systems in any given galaxy are ex- 
pected to exhibit a range of fo values. For example, disrup- 
tion rates of GCs in dense high-redshift environments were 
likely higher than in current quiescent galaxies (Kruijssen et 
al. 2012) due to the enhanced strength and rate of tidal per- 
turbations. Likewise, GCs on eccentric orbits with smaller 
perigalactic distances R per i should have smaller fo values than 
otherwise (initially) similar GCs with larger values of R ve „. 
Finally, the rate of mass loss from GCs by evaporation, /z ev , 
effectively scales with the mean GC half-mass density ph as 
Mev oc p" 5 for tidally limited GCs (e.g., McLaughlin & Fall 
2008; Gieles et al. 2011; Goudfrooij 2012). Among GCs for 
which two-body relaxation has been the dominant mass loss 
mechanism, high-ph GCs thus feature smaller fo values than 
low-ph GCs at a given M c \- 

Figure 7 shows the present-day MF slopes oimf in the stel- 
lar mass range 0.3-0.8 M Q according to the K09 models 
for surviving GCs as a function of A4 c i- The two panels (a) 
and (b) show the results for two different slopes of the stellar 

6 to is defined by ta s = to M"l\ where ?di s is the cluster disruption time and 
7 sets the mass dependence of cluster disruption. See K09 for details. 



IMF: olimf = -3.00 and -2.35, respectively. We also indicate 
"weighted mean" values of cxmf for a GC population with 
M c \ > 2.5 x 10 5 M Q at an age of 12 Gyr that was formed 
with a Schechter (1976) initial cluster mass function and a 
maximum initial mass of 10 7 M Q . For reference, a power-law 
fit to the Kroupa IMF in the mass range 0.3-0.8 M Q yields 
oluf = — 1.7 (cf. panels (d) and (h) of Fig. 1). 




5 5.5 6 6.5 

Log (M c \ /M ) at age = 1 2 Gyr 



FIG. 7. — Panel (a): MF slope ot M F in the range 0.3 < M-^/Mq < 0.8 
versus GC mass at an age of 12 Gyr for ce/Mf = -3.00 according to solar 
metallicity models of K09 that produce surviving GCs with M c \ > 2.5 X 
10 s Mq at an age of 12 Gyr. Solid lines represent models with King (1966) 
parameter Wo = 7, while dashed lines do so for Wo = 5. From top to bottom, 
the lines represent models with to = 3.0, 1.0, 0.6, and 0.3 Myr (the latter 
only for Wq = 7). Circles indicate weighted mean values for the abscissa and 
ordinate for each model, assuming a Schechter (1976) initial cluster mass 
function to assign weights and cluster masses 2.5 X 10 5 < M c i/M e < 10 7 
at an age of 1 2 Gyr. The magenta dash-dotted line (along with an error bar) 
indicates the weighted average olmf of all circles (using inverse variances as 
weights), while the green dotted line indicates the Kroupa IMF. Panel (b): 
Similar to panel (a), but now for a Salpeter (1955) IMF (oimf = -2.35). 

Focusing on the "weighted mean" values of a Mf for the 
K09 models shown in Figure 7, we find overall mean values 
omf = -1 .84 ± 0.26 for a IMF = -3 .00 and = -1 .80 ± 0. 1 1 
for climf = -2.35, indicating that GCs with a steeper IMF ex- 
perience a stronger evolution of the stellar MF. Note that both 
of these values are consistent with a Kroupa MF and hence 
consistent with the observed mean colors of red GCs shown 
in Figure 1 . It is important to realize that the stronger flat- 
tening of olmf over time in massive GCs with steeper IMFs is 
mainly caused by the weaker effect of retained massive stel- 
lar remnants on the escape rate of massive stars in GCs with 
steeper IMFs relative to that in GCs with flatter IMFs. As 
described in detail in K09, retained massive stellar remnants 
boost the escape rate of massive stars relative to that of low- 
mass stars. For a Kroupa IMF, this effect limits the flatten- 
ing rate of olmf over time relative to GCs with a lower frac- 
tion of massive stellar remnants (see K09). However, GCs 
with steeper IMFs feature a smaller mass fraction taken up by 
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massive remnants. This notably reduces the 'dampening ef- 
fect' of retained massive remnants on the evolution of czmf- 
As a result, the present-day stellar MFs in GCs are consis- 
tent with any IMF slope a < -1.7. This degeneracy implies 
that the IMF cannot be uniquely constrained in this range for 
individual GCs. However, the large spread of observed MF 
slopes among Galactic GCs (de Marchi et al. 2007; Paust et 
al. 2010) seems to be more consistent with a bottom-heavy 
IMF (compare both panels in Figure 7). 

We conclude that the difference in g—z orB-7 color between 
giant elliptical galaxies and their "average" massive red GCs 
shown in Figures 1 and 2 is consistent with long-term dynam- 
ical evolution of star clusters that were formed with a steeper 
IMF than those derived from star counts in the solar neighbor- 
hood. Among the three mechanisms we think could conceiv- 
ably cause this color difference, this is the only one that seems 
consistent with the data. Hence we suggest that this color 
difference constitutes evidence for a "bottom-heavy" IMF in 
luminous elliptical galaxies. 

6. DISCUSSION 

Our results add to a developing consensus that luminous el- 
liptical galaxies have "heavy" IMFs, with a greater contribu- 
tion from low-mass dwarf stars than IMFs determined from 
star count studies in the Milky Way. Cenarro et al. (2003) 
studied the Call triplet near 8600 A, which is stronger in 
RGB stars than in MS stars, in a sample of elliptical galax- 
ies. They found that the strength of the Ca II triplet decreases 
with increasing velocity dispersion among elliptical galax- 
ies. This is hard to explain in terms of an Ca abundance 
trend since Ca is an a element, and [a/Fe] as derived from 
Mg and Fe features near 5200 A does not show this trend. 
Hence, they suggested that this trend might reflect an increas- 
ing IMF slope with increasing galaxy mass, as indicated by 
SSP models of Vazdekis et al. (2003). Similarly, van Dokkum 
& Conroy (2010, 2011) found strong equivalent widths of 
the Nal AA8183,8195 doublet and the molecular FeH band 
at 9916 A in near-IR spectra of giant elliptical galaxies. As 
both of those indices are stronger in red dwarf stars than 
in red giants, their strongly favored explanation is the pres- 
ence of a significantly larger number of late M dwarf stars 
than that implied by a Kroupa IMF, even though there is a 
slight degeneracy with [Na/Fe] ratio (see also Conroy & van 
Dokkum 2012; Smith et al. 2012). Finally, Cappellari et al. 
(2012) studied stellar M/L ratios [(M/L) stMS ] derived from 
integral-field spectroscopy of 260 early-type galaxies. They 
found that the ratio of {M/L) Am& by the M/L predicted by 
SSP models that use a Salpeter IMF [(M/L)s d ip] varies sys- 
tematically with galaxy velocity dispersion: galaxies with the 

7 If the present-day MF of those GCs is flatter than the Kroupa or Chabrier 
IMFs, as might be expected given the strong tidal field at their location deep 
in the Galaxy's potential well, their optical colors would be bluer than if they 
had Kroupa or Chabrier MF's (cf. Figures 1 and 2). Thus, absolute ages 
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largest values of (M/L) st . drfi /(M/L)s. d i p have the highest ve- 
locity dispersions (see also Tortora et al. 2012). However, 
dynamical modeling such as that done by Cappellari et al. 
(2012) cannot unambiguously constrain whether the highest 
values of (A / J/L) st ars/(-A / J/L)saip are due to a large population 
of low-mass stars (i.e., a bottom-heavy IMF) or large number 
of massive stellar remnants (i.e., a top-heavy IMF). The same 
degeneracy is present for mass estimates from gravitational 
lensing (e.g., Treu et al. 2010). In contrast, our results indi- 
cate that the bottom-heavy IMF is the more likely solution. 

Finally, we review the dependence of our results on SSP 
model ingredients. At old ages (> 10 Gyr) and high metal- 
licities ([Z/H] > -0.5), SSP models still suffer from several 
limitations. The likely most significant internal uncertainty 
is related to lifetimes of AGB stars (see, e.g., Melbourne et 
al. 2012). In terms of external uncertainties, we note that 
integrated colors and spectral indices of several SSP models 
are calibrated to observations of the two Galactic bulge GCs 
NGC 6528 and NGC 6553 that have Z ~ Z (e.g., Thomas et 
al. 2003; Maraston 2005). However, there are no published 
MFs for those two GCs available to our knowledge, yielding 
an uncertainty in the absolute calibration 7 . Notwithstanding 
such limitations, our main result is based on g—z or B-I color 
differences between the average red massive GC and their par- 
ent galaxies and the absence of such differences in H/3 and 
[MgFe]' indices. As the interpretation of relative colors or 
line strengths should be largely model-independent, our re- 
sult strongly suggests a bottom-heavy IMF in massive ellipti- 
cal galaxies. 
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and/or [Z/H] values of old metal-rich populations derived from integrated 
colors would currently be slightly overestimated. 
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APPENDIX 

BUILDING SYNTHETIC INTEGRATED COLORS OF METAL-RICH GLOBULAR CLUSTERS CONTAINING TWO STELLAR 

GENERATIONS 

In this appendix, we determine colors of various stellar types from synthetic spectra with chemical compositions typical of first— 
and second-generation stars found in massive metal-rich GCs. We assume [Z/H] = 0.0 and explore two choices for the abundances 
of He, C, N, O, and Na. To simulate first-generation (FG) stars, we choose the primordial He abundance (Y = 0.235 + 1.5 Z), we 
follow Cannon et al. (1998) by choosing [C/Fe] = 0.06 and [N/Fe] = 0.20, and we adopt [O/Fe] = 0.40 and [Na/Fe] = 0.00 from 
Carretta et al. (2009). To simulate second-generation (SG) stars, we choose [C/Fe] = -0.15, [N/Fe] = 1.05, [O/Fe] = -0.10, and 
[Na/Fe] = 0.60 (see Carretta et al. 2009). For the He abundance of SG stars, we consider the case of NGC 6441, a metal-rich GC 
in the bulge of our Galaxy ([Fe/H] = -0.59, Zinn & West 1984). NGC 6441 (as well as the similar cluster NGC 6388) is unusual 
among metal-rich GCs in our Galaxy in that it shows a horizontal branch (HB) that extends blueward of the RR Lyrae instability 
strip in the CMD whereas all other metal-rich GCs have purely red HBs (e.g., Rich et al. 1997) 8 . The blueward extension of the 
HB of NGC 6441 is widely thought to be due to enhanced He in a fraction of its constituent stars, and we adopt the value Y = 0.33 
from Caloi & D'Antona (2007). As blue HBs are very unusual among massive metal-rich GCs in our Galaxy, this value of Y 
should probably be considered an upper limit for SG stars in the context of the current paper. 

To produce synthetic integrated colors of stellar populations, we first build two sets of synthetic spectra for six stellar types that 
span almost the full range of luminosities L and temperatures T e ff encompassed by the isochrones (see Figure 8). For this purpose 

8 Note however that the impact of the blue extension of the HB seen in NGC 6388 and NGC 6441 to the overall optical color of their HB is negligible. In 
terms of A(V-7), defined by Dotter et al. (2010) as the median color difference between the HB and the RGB at the luminosity of the HB, the averaged HB of 
NGC 6388 and NGC 6441 is 0.003 ± 0.010 mag redder than the average HB of the 8 other Galactic GCs with Age > 12 Gyr and -0.70 < [Fe/H] < -0.50 in that 
study, i.e., Lynga 7, NGC 104, NGC 5927, NGC 6304, NGC 6496, NGC 6624, NGC 6637, and NGC 6838. 
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TABLE 2 

Stellar parameters of the model atmospheres. 



r eff (fg) 


r eff (so) 




Vturb 


(i) 


(2) 


(3) 


(4) 


3602 


3665 


4.79 


2.0 


4526 


4795 


4.59 


2.0 


5630 


5700 


4.20 


2.0 


4835 


4941 


3.48 


2.0 


4237 


4311 


1.96 


2.0 


3537 


3608 


0.70 


2.0 



NOTE. — Column (1): Effective temperature for first-generation stars in K. Column (2); Effective temperature for second-generation stars in K. Column (3): 
logarithm of surface gravity in cm s -1 . Column (4): Turbulent velocity in km s" 1 . 

we select Dartmouth isochrones 9 (Dotter et al. 2007) with [Z/H] = 0.0, [a/Fe] = +0.4, and the two values of Y mentioned in the 
previous paragraph. [Z/H] values for alpha-enhanced populations are evaluated using [Z/H] = [Fe/H] + 0.929 [a/Fe] (Trager et al. 
2000). The parameters of the model atmospheres are listed in Table 2. These temperatures, gravities, and chemical abundances 
are then used to calculate model atmospheres and synthetic spectra using the codes ATLAS 12 and SYNTHE (Kurucz 2005; 
Castelli 2005; Sbordone et al. 2007), respectively. ATLAS12 allows one to use arbitrary chemical compositions. In doing so, we 
use model atmospheres from the [Fe/H] = 0.0, [a/Fe] = +0.4 model grid of Castelli & Kurucz (2004) 10 as reference and adjust 
Tea, g, [Fe/H], and the individual element abundances in the ATLAS 12 runs. The resulting synthetic spectra are then redshifted 
according to the heliocentric radial velocity of the Virgo galaxy cluster (1035 kms" 1 , Mould et al. 2000) and integrated over 
passbands of Johnson/Cousins (Uj, Bj, Vj, R c , lc, Bessell 1990) and Sloan ugriz 11 to produce synthetic magnitudes and colors. 




log (Te„) 

FIG. 8. — Theoretical isochrones from the lower MS to the upper RGB for the two different chemical compositions associated with two stellar generations 
within massive GCs. The solid line corresponds to the first-generation population (with primordial He abundance), whereas the dashed line corresponds to the 
second-generation population (with enhanced He abundance Y = 0.33). The filled and open squares indicate the stellar types for which we calculated model 
atmospheres and synthetic spectra (cf. Table 2). See text for more details. 

Three examples of the synthetic spectra for FG and SG stars are shown in Figs. 9-11. Labels in each plot indicate prominent 
absorption features that change significantly in strength between the two chemical compositions. These spectra are described 
briefly below. 

Figure 9 compares RGB spectra of the FG and SG mixtures. The higher [N/Fe] of the SG mixture causes stronger molecular 
bands of CN. This is especially clear at A > 7000 A. The main reason why this is not as obvious for the CN bands near 3590, 
3883, and 4150 A and the NH band near 3360 A is that the enhanced absorption in the blue region for the SG mixture is partially 
compensated by its enhanced He abundance which elevates the stellar continuum in the blue. The increased opacity at A < 4500 
A for the SG mixture also causes a somewhat elevated stellar continuum at longer wavelengths, which likely causes the slightly 
higher flux of the SG mixture in the 4500-6200 A region. Finally, the lower O abundance of the SG mixture causes a weaker 
OH absorption band between 3050 and 3200 A which just falls within the Uj and u bands for the radial velocities of the galaxies 
in our sample. 

In the spectra of TO stars (Figure 10), the higher temperature causes generally weaker molecular bands, rendering generally 

9 see http://stellar.dartmouth.edu/ models/index.html 

10 see http://wwwuser.oat.ts.astro.it/castelli/grids.html 

11 see http://www.sdss.Org/drl/instraments/imager/index.html#filters 
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FIG. 9. — Top panel: Synthetic spectra for RGB models with log g = 1.96 for the FG (black) and SG (red) mixtures. The spectra were convolved with a 
Gaussian with a = 12 A to improve readability. Relevant molecular bands are labelled. Bottom panel: Transmission curves for Johnson/Cousins Uj, Bj, Vj, Rc, 
and If filters (solid curves, left to right) and SDSS ugriz filters (dashed curves). See text for details. 




small differences between the FG and SG spectra. Only the NH and CN bands remain slightly stronger in the SG mixture. 

The differences between the FG and SG spectra of cool MS stars (Figure 1 1) are generally similar to those seen for the RGB 
spectra, except that the MS stars show weaker red CN bands and a slightly stronger OH band in the 3000-3200 A region. 

Differences in stellar absolute magnitude AM = Msg -Mfg are then evaluated as a function of log g for all filter passbands 
considered here. The results are illustrated in Figure 12 which shows AM as a function of log g. The largest differences show 
up in the Uj and u passbands, due to their sensitivity to relevant parameters: (i) T e ff and hence Y; (ii) N abundance through 
the NH band near 3360 A as well as the CN bands near 3590 and 3883 A; and (Hi) O abundance through the OH band at the 
short-wavelength edge of the U j and u filters. The weaker OH absorption for the SG population is the main reason why AM for 
the Uj and u passbands is negative at the top of the RGB (i.e., lowest logg). AM then increases significantly with increasing log 
g to become positive at the TO (due mainly to stronger NH and CN absorption in the SG population), and then decreases again 
to the bottom of the MS (i.e., highest log g). The other passbands (longward of ~ 4000 A) show a similar dependence of AM on 
log g, although the absolute AM values at A > 4000 A are much smaller than for U j and u at a given log g. 

The dependence of AM between FG and SG populations as a function of wavelength found here are qualitatively similar to 
those of Sbordone et al. (2011) who did a similar study at low metallicity ([Fe/H] = -1.62). However, we find that AM values 
are different for metal-rich populations in a quantitative sense, sometimes by a significant amount. 

Isochrone tables for the SG population are then calculated by adding AM to the absolute magnitudes for each filter passband 
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FIG. 12. — Differences between stellar absolute magnitudes of the FGand SG stars considered in this paper as a function of log g. Upper panel: Johnson/Cousins 
UBVRI filters (see legend for symbols). Lower panel: Sloan ugriz filters. Dashed vertical lines delineate the RGB, TO, and MS regions on the CMD (see labels 
near the bottom of the upper panel). See Appendix A for more details. 

in the (FG) isochrone tables (i.e., Msg = M?g + AM), using linear interpolation in log g. After rebinning the resulting isochrone 
tables to a uniform bin size in stellar mass using linear interpolation, weighted integrated luminosities of the population were 
derived for the relevant passbands by weighting individual stellar luminosities in the isochrones using a "standard" Kroupa 
(2001) IMF. Resulting integrated-light magnitude offsets (in the sense SG - FG) are listed for each filter passband in Table 3. 
Interestingly, the large differences found in the Uj and u passbands between FG and SG populations at a given stellar type largely 
cancel out in integrated light (i.e., when integrated over the stellar mass function). Integrated-light magnitude offsets are found 
to stay within 0.03 mag in an absolute sense for any of the filter passbands considered here. Finally, and in the context of the 
subject of this paper, we note that the resulting offsets between SG and FG populations in the B—I and g-z colors are -0.041 
and -0.054 mag (in the sense SG-FG), respectively. Since recent multi-object spectroscopy studies of massive Galactic GCs 
typically show an approximate 50-50% split between FG and SG stars (e.g., Gratton et al. 2012), we conclude that the predicted 
overall effect of the presence of multiple stellar generations on the integrated B—I and g-z colors of massive metal-rich GCs is 
offsets of -0.021 and -0.027 mag, respectively. 
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TABLE 3 

INTEGRATED-LIGHT MAGNITUDE OFFSETS BETWEEN THE SG AND FG POPULATIONS. 



AU, 


ABj 


AVj 


AR C 


AI C 


+0.002 


-0.018 


-0.027 


-0.002 


+0.023 


Aw 


Ag 


Ar 


a; 


Az 


+0.022 


-0.031 


-0.007 


+0.018 


+0.023 



NOTE. — Magnitude offsets are in the sense SG - FG. 



